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ABSTRACT 

We introduce methods to quantify the X-ray morphologies of supernova rem- 
nants observed with the Chandra X-ray Telescope. These include a power-ratio 
technique to measure morphological asymmetries, correlation-length analysis to 
probe chemical segregation and distribution, and wavelet-transform analysis to 
quantify X-ray substructure. We demonstrate the utility and accuracy of these 
techniques on relevant synthetic data. Additionally, we show the methods' capa- 
bilities by applying them to the 55-ks Chandra ACIS observation of the galactic 
supernova remnant W49B. We analyze the images of prominent emission lines 
in W49B and use the results to discern physical properties. We find that the 
iron morphology is very distinct from the other elements: it is statistically more 
asymmetric, more segregated, and has 25% larger emitting substructures than 
the lighter ions. Comparatively, the silicon, sulfur, argon, and calcium are well- 
mixed, more isotropic, and have smaller, equally-sized emitting substructures. 
Based on fits of XMM-Newton spectra in regions identified as iron rich and iron 
poor, we determine that the iron in W49B must have been anisotropically ejected. 
We measure the abundance ratios in many regions, and we find that large, lo- 
cal variations are persistent throughout the remnant. We compare the mean, 
global abundance ratios to those predicted by spherical and bipolar core-collapse 
explosions; the results are consistent with a bipolar origin from a ~25 Mq pro- 
genitor. We calculate the filling factor of iron from the volume of its emitting 
substructures, enabling more precise mass estimates than previous studies. Over- 
all, this work is a first step toward rigorously describing the physical properties 
of supernova remnants for comparison within and between sources. 
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1. Introduction 

Supernova remnants (SNRs) are a diverse class of objects that play an essential role 
in the dynamics of the inte rstellar medium (ISM) an d that produce and distribute most of 



the metals in the Universe (IFukugita fc Peebles! l2004l ). Since its launch, t he Chandra X-ray 



Telescope has observed over one- hundred SNRs in the Milky Way galaxy (lGreedl2006l ). The 
sub-arcsecond spatial resolution and spectro-imaging capabilities of Chandra have facilitated 
detailed studies of the metal-rich ejecta from supernova explosions as well as the interactions 
with their environments as they expand (see Weisskopf & Hughes 2006 for a review). Chandra 
has revealed that X-ray spectral and spatial morphologies of SNRs are extremely complex. 
However, no systematic, quantitative methods for comparison within and between SNRs 
have been applied to develop a global picture of these varied sources. 

Rigorous mathematical techniques are necessary to maximally utilize these Chandra 
datasets. Toward this end, we apply well-estalished tools to characterize the X-ray mor- 
phologies of SNRs. We present three methods especially suited for application to Chandra 
Advanced CCD Imaging Spectrometer (ACIS) observations: a multipole power-ratio tech- 
nique to probe morphological symmetry, correlation-length analysis (an adaptation of two- 
point correlation) to measure chemical segregation and distribution, and wavelet-transform 
analysis to quantify substructure. These methods have been used previously in other astro- 
nomical contexts, and we have extended and further developed them for SNR science. In 
§2, we describe the methods in detail and demonstrate their application to synthetic data of 
extended sources. In §3, we apply these techniques to the 55-ks Chandra ACIS observation 
of W49B, a SNR with a particularly complex spatial distribution. In §4, we discuss how the 
results can be used to discern between the heating and explosion properties of the source. 
We present our concluding remarks in §5. 



1.1. W49B Background 

W49B (G43.3— 0.2) is a principal example of a galactic SNR with complex X-ray mor- 
phology. It is the brightest ejecta-dominated SNR in X-rays {Lx ~ 10^^ erg s~^; Immler 
& Kuntz 2005), and it has one of the highest surface brightnesses at 1 GHz of all galactic 
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sources (IMoffett &: Reynoldslll994l ). Radio maps of the source show W49B has a prominent 
shell structure ~4' in diameter, yet it has a center- filled morph ology in the X-rays. HI ab - 
sorption measurements indicate the remnant is 8±2 kpc away (IMoffett &: Reynolds! 1 19941 ). 
and due to its distance through the Galactic plane, no optical emission is detected from the 
source. A recent Spitzer study of SNRs in the Galactic plane observed filament ary structure 



of inf rared line emission in W49B from both shocked ionic and molecular gas (IReach et al. 
2006h . 



W49B was first detected in X-rays by the Einstein Observatory ( iPye et al.lll984l ). An 
EXOSAT observation soon after revealed intense iron line emission from the source. ASCA 



spect ra confirmed the high abundances of several metals in W49B ( Pye et al.lll984l : lHwang et al, 
20001 ). and the emission was described best by two distinct thermal components of optically- 
thin plasma. Given the remnant's relatively young age (~ 1000 years: Pye et al. 1984; Smith 
et al. 1985), the line emission was attributed primarily to the ejecta from the progenitor 
star. Recent hi gh-resolution X-ray images (see Figure [1]) and spectra from Chandra and 
XMM-Newton (IMiceli et al.l l2006l : iKeohane et al.l 120071 ) have shown that the metals have 
remarkable spatial distributions. The emission is elongated in a centrally-bright bar that is 
enriched in iron and nickel and has two plumes at its Eastern and Western edges that are 
perpendicular to the axis of this bar. The Eastern part of the remnant terminates sharply, 
while the Western portion of W49B is remarkably diffuse. Additionally, the Western regions 
have weak iron ernission , whereas the lower-Z elements have significant emission in these 
areas JlVIiceh et al.lbood ). 



Two potential scenarios may account for the complex morphology of W49B (Miceli et al. 
2006, 2008; Keohane et al. 2007). The high-metallicity and barrel-shaped morphology could 
be signatures of a bipolar explosion (possibly a hypernova) of a massive star. Conversely, 
a complex environment around a typical, symmetric supernova explosion could give rise to 
inhomogeneous heating to X-ray emitting temperatures. The latter is an appealing model 
as it is a commonly observed phenomenon (e.g., SNR G292. 8-1-1. 8: Park et al. 2004). 



However, this scenario would require unrealistic ejecta masses (IKeohane et al.l 120071 ). The 
bipolar origin of W49B is probl ematic as well because the predicted rate of such events is 
only 1 per 10^ years per galaxy (jPodsiadlowski et al.ll2004l ). Thus, it is uncertain presently 
which scenario is the correct interpretation of W49B's exotic shape. 

In order to break the degeneracy of environment versus explosion effects, quantitative 
methods can be applied to the existing data. The Chandra images of W49B contain a wealth 
of physical information that can be extracted using appropriate mathematical anaylses. Mea- 
surement of elemental segregation and mixing are important clues to explosion and evolu- 
tion histories. Additionally, a rigorous description of X-ray substructure enables much more 
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reliable estimates of elemental abundances and temperature variations than any previous 
work. Together, these calculations can distinguish the correct interpretation of asymmetries 
(anisotropic ejection o r heating inhora ogeneities), one of the biggest outstanding questions 
in SNR studies today (ICanizaresI 120041 ) . 

The anomalous spatial distribution and substructure of W49B make it an ideal source 
to test our methods. We aim to define quantitatively the spatial and spectral features of 
W49B as a means to extract physical insights about this enigmatic source. 



2. Quantitative Methods for Imaging Analyses 

In this section, we describe the techniques that can be used to quantify X-ray morphol- 
ogy of extended sources. We present their mathematical formalism as well as examples of 
their application to synthetic data. 



2.1. Power- Ratio Method 

The power-ratio method (PRM) was developed initially to quantify the X-ray morphol- 
ogy of galaxy clusters observed with ROSAT as a probe of their dynamical states (Buote 
& Tsai, 1995, 1996). The technique was extended to Cha ndra observations to measure the 



evolution of cluster morphology as a function of redshift (j Jeltema et al.l l2005l ) . The PRM 
measures sensitively asymmetries in an image via calculation of the multipole moments of 
the X-ray surface brightness in a circular aperture. We present an overview of this technique 
below, and we refer the reader to the above papers for the rigorous details and examples of 
its application to synthetic data and X-ray observations. 

The PRM is derived similarly to the multipole expansion of the two-dimensional gravi- 
tational potential interior to radius R: 

^'(i?, 0) = -2Gao\n f — j - 2(7 ^ — ^ (a^cosm^ + fe^sinm^) . (1) 

^ ^ m=l 

where the moments am and hm are 

am{R) = I S(f') W"cosm0'rf"x 
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x' — {R', (f)'), and E is the surface mass density. For our imaging analyses, the X-ray surface 
brightness replaces surface mass density in the power ratio calculation. 

The powers of the multipole expansion are obtained by integrating the magnitude of 
"ifm (the mth term in the multipole expansion of the potential) over a circle of radius R, 

Pm{R) = ^ r *-(^' '^)*-(^' '^)^'^- (2) 

Ignoring the factor of 2G, this equation reduces to 



Po = [aoln(i?)]' (3) 

^ + ■ 

The moments and hm (and consequently, the powers Pm) are sensitive to asymmetry 
in the X-ray surface brightness distribution, and the higher-order terms measure asymmetries 
at successively smaller scales relative to the remnant size. To normalize with respect to flux, 
we divide the powers by Pq to form the power ratios, Pm/Po- Pi approaches zero when the 
aperture is centered on the centroid of an image, so we have set the origin in all analyses to 
the centroid of each individual image. In this case, morphological information is given by 
the higher-order terms P2/PQ, P3/P0: and Pi/Po- The quadrupole power P2 is sensitive to 
the ellipticity of a source. The moment P3 measures deviations from mirror symmetry (such 
as triangular morphology), and it is not affected by mirror-symmetrical features. The power 
P4 is correlated with P2 but is more sensitive to smaller-scale structure than P2. Thus, the 
power ratios provide complementary measures of X-ray morphology in extended sources. 



2.2. Correlation-Length Analysis 

Two-point correlation has been applied to measure galaxy clustering (e.g., Davis & 
Peebles 1983) as well as to map the anisotropics of the Cosmic Microwave Background 
(e.g., Kogut et al. 1996). We have adapted this method to characterize spatial distribution 
and segregation between elements in an X-ray source, and we call our modified technique 
correlation- length analysis (CLA). 

CLA measures essentially whether two images have statistically similar spatial distribu- 
tions. We define the correlation length of each pixel in one image as the minimum distance 
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(>0 pixels) to another pixel of equal intensity in a second image. The two images are nor- 
malized to have the same maximum intensity, and pixel values are rounded to give integer 
counts throughout the images. We calculate the correlation length for each pixel, and we 
create maps of these values to identify regions with and without particular emission features. 

Additionally, we determine whether two images have similar overall spatial distributions 
by calculating the difference between the cumulative distributions of their correlation lengths 
to the cumulative distributions of each of their auto-correlations. If the two images have 
distinct spatial distributions, the correlation lengths between them will be disparate from 
those of their auto-correlations. Thus, the difference in the trends of correlation length 
values provides a measurement of the similarity of two images. We perform auto-correlations 
on mock images produced via a Monte Carlo approach to estimate the uncertainty of the 
calculation. 

We test the utility of CLA by applying the method to three groups of ten synthetic im- 
ages (Figure [2] has examples and a mathematical description of the datasets). Groups 1 and 2 
have similar distributions and represent auto-correlations, while Group 3 reflects comparison 
of disjoint distributions. We performed CLA on every combination (ninety different pairings) 
of the ten images within each group and produced the three mean correlation-length CDFs 
and their 1-cr ranges shown in Figure [31 

To determine whether the size of substructure influences CLA results. Groups 1 and 

2 have clumps at fixed location but with scales selected from Gaussian distributions of 
different variance: the clump-size distribution of Group 2 has 56% larger variance than that 
of Group 1. Comparing their resulting CDFs, we find that the confidence ranges overlap 
for these datasets, with an average separation of 0.5cr. Thus, images in Groups 1 and 2 
have statistically similar spatial distributions, indicating that clump-size variation does not 
hinder CLA from distinguishing homogeneous populations. 

To ascertain how clump position affects CLA results, clump locations of Group 3 
datasets were selected randomly while the clump sizes did not vary much (with 11% of 
the Group 1 variance). In this case, many pixels have large correlation lengths, so the Group 

3 CDF is very disparate from those of Groups 1 and 2 with average separations of 3. la 
and 2.7(7, respectively. Thus, Group 3 has a statistically distinct spatial distribution from 
Groups 1 and 2, and we conclude that CLA is a good measure of whether substructures 
have similar locations. Therefore, CLA is well-suited to compare morphologies within and 
between extended sources. 



- 7- 



2.3. Wavelet-Transform Analysis 



Wavelet-transform analysis (WTA) was applied successf ully to ROSAT and E instein 
data to extract small-scale X-ray structure of galaxy clusters (iGrebenev et al.lll995l ). With 
the superior spatial resolution of Chandra, this technique can be extended to characterize 
precisely the size and distribution of emitting regions in X-ray sources. A wavelet basis 
enables scale as well as spatial localization, thus WTA can identify substructure, its extent, 
and its position. 



2.3.1. Mathematical Formalism 

A wavelet-transformed image is a decomposed image of a signal's intensity (from herein, 
power) measured at the scale of a filter size. Mathematically, a wavelet transform w is the 
correlation of a signal s{x,y) in an image with the analyzing wavelet function g{x,y): 

w{x,y,a) = s{x,y) ^-g(-,-), (5) 

a \a a/ 

where a is the scale (or width) of the wavelet transform. For astronomical images, a ra- 
dial Mexican-hat function (the normalized second derivative of a Gaussian function) is the 
optimal wavelet to use because of its similar shape to Gaussian signals an d it removes fiat 



features (like diffuse emission) since it has a zero mean (jSlezak et al.lll990l ). In the analysis 



that follows, we use a radial Mexican-hat wavelet qi-,-) of the form 



' a a' \ 



Wavelet-transformed images are produced by calculating w{x, y, a) for each pixel (m, n) 
in a raw image: 

/ \ 1 \ ^ f-^i -^m yj yn\ 

w{m,n,a) = -} CijQl , , (7) 

where Cij is the number of counts in the {i,j) pixel. If we assume the number of counts in 
a given pixel is independent of the counts in nearby pixels, we can use Poisson statistics to 
calculate the variance q{m, n, a) of w{m, n, a): 



q{m, n,a) = ^ijg^ (— — 



2 j Xi Xjji yi 

5 
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where dij is the variance of Cij in the pixel of a raw image. 

Essentially, w measures the summed intensity enclosed by the area of the Mexican hat. 
Thus, the size of an individual source can be characterized by the scale where the convolution 
of the wavelet and a signal reaches a maximum. To exemplify this point mathematically, we 
assume that a source has an isotropic Gaussian shape with signal s{x,y,a) given by 



s{x,y,a) = ^^^exp — , (9) 



where / is the intensity of the source located at (xq, i/o). Using Eq. 5 and Eq. 6, the wavelet 
coefficient w{xo,yo,ci) at this position is 

21 ( 

w{xo,yo,a) = — 1 + ^ • (10) 
If we divide Eq. 10 by a, we find that 

W 2// ^2 2 

- = —! + — • (11) 

Eq. 11 has a global maximum at Omax = o". Consequently, by identifying the peak wja value 
at the central pixel of an emitting region, we measure its characteristic size. 



2.3.2. Application to Synthetic Data 

To demonstrate the utility of WTA for astronomical images, we apply the method to 
relevant synthetic datasets. We begin by using WTA to extract the size of a single two- 
dimensional Gaussian function (from herein, a clump) with a width cr. Figure UK shows this 
synthetic data, and Figure Hb gives the plot oiw/a versus a (intensity per unit cLrG3i clS 3i 
function of scale; from herein, WTA plot) for the central pixel of the clump. Larger w/a 
values correspond to more power emitted at the given scales, so local and absolute maxima 
indicate the size of emitting regions. Indeed, the peak Omax of the WTA plot occurs at the 
exact value of a. Therefore, the method can measure accurately the size of an individual 
emitting region or substructure. 

By applying WTA to all the pixels in an image and summing their results for each 
scale, the method can be utilized to calculate the intensity profile of an entire extended 
source. As an example, we have applied WTA to all pixels in three synthetic datasets with 
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randomly-distributed clumps of equal width (Figure [5|, left). The number of clumps in these 
images are varied to demonstrate the effects of increased emitting surface area, or filling 
factor. Figure O (right; solid lines) displays the averaged WTA plot {{w)/a versus a) for 
the three datasets. The curves have similar shape as Figure IH and increased {w)/a values 
correspond to more power emitted at a particular scale over an entire image. As the number 
of clumps increases and they overlap, larger agglomerates are formed. This effect increases 
the surface area of contiguous emitting regions, and the amount of power measured at large 
scales to rise. Thus, increased numbers of clumps shift the peak amax to higher values, and 
the turnover in {w)/a becomes less prominent. 

In images with many clumps, we can measure the scale of isolated structures using the 
WTA plot for their central pixels. However, larger filling factors cause individual clumps 
to agglomerate, and the number of segregated structures decrease as a result. In Figure [5] 
(right; dashed lines), we plot the averaged WTA plot for the individually identified clumps 
in the three synthetic datasets of randomly-distributed clumps of equal width (Figure [5l 
left). We define the centroid of individual clumps as the locations with local maxima in 
w/a separated at least three pixels from another local maximum. In all cases, the algorithm 
successfully reproduces the scale of the isolated structures. Generally, our ability to identify 
individual clumps decreases with greater emitting surface area (see Figure [5]). Consequently, 
WTA probes an extended source's filling factor as well as the size of its individual emitting 
regions. 

Real astronomical images have noise that will naturally influence these WTA results. To 
investigate the effects of noise, first we apply our techniques to synthetic data of pure noise 
(Figure [6l left). The noise intensity in each pixel was selected from a Poisson distribution. 
Figure [6] (right) shows the resulting averaged WTA plot; the absolute maximum of the WTA 
curve is at a = 1 pixel. Since the noise is distributed on a pixel-by-pixel basis, it is natural that 
the most power is emitted at the single-pixel scale. We then add the pure noise (Figure El 
left) to synthetic data of randomly-distributed clumps to simulate realistic images of an 
extended astronomical source. Clump sizes were selected from Gaussian distributions, and 
images were produced using different variances in the clump-size distribution to investigate 
the effects of substructure on WTA results. We calculate the averaged WTA plot for the 
datasets; two example images and their results are shown in Figure [71 All curves have 
global maxima at Omax = 1 pixel since the noise is dominant. However, we can reduce its 
contribution by ignoring all pixels with global maxima at the smallest scales. Using this 
approach, the global maxima shift toward the scale of mean clump size since the most power 
resides there when noise is excluded. We determine that the variance in clump size did 
not alter the WTA results either: the intensity profile is similar regardless of how divergent 
individual clumps are from the average scale. 
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Instrumental effects may alter WTA results: in particular, an image's resolution is lim- 
ited by the point-spread function (PSF) of the detector. To investigate how an instrument's 
PSF changes the scale where most power is emitted, we convolved the data in Figure [7^ with 
Gaussian functions and then add noise. The resulting four images are shown in Figure [HI 
along with their averaged WTA plot. Generally, the convolution did not influence the inten- 
sity profile measurement. We find that the results are unaffected for all PSFs 60% or less of 
their scale, and the maxima in the WTA plot increase with PSFs larger than this fraction. 

Based on the results from the synthetic data presented above, we conclude that WTA is 
a useful technique to extract the individual sizes of emitting regions and to measure a source's 
filling factor, an important geometric parameter. Noise increases the power extracted at the 
smallest scales, but we can limit this contribution by ignoring pixels with global maxima 
of a single pixel. Generally, large variation in clump sizes does not reduce the accuracy in 
determining their intensity profiles, and WTA measures the various emitting regions' scales 
necessary to produce a given surface brightness. Thus, WTA is a useful tool to probe the 
substructure of extended sources, such as supernova remnants. 



2.4. Uncertainties and Limitations of the Methods 

Although the methods are successful on these synthetic data, the techniques have un- 
certainties and limitations that should be considered. The power-ratio method necessitates 
an accurate determination of an image's centroid; this calculation requ ires a source to have a 



minimum flux and minimum signal-to- noise ratio (S/N). As shown by lJeltema et al.l (120051 ). 
circular apertures with greater than 500 net counts and S/N greater than 3.0 are sufficient 
to locate the centroid for unbinned Chandra data. Furthermore, to obtain accurate net 
moments of extended sources, a few criteria should be satisfied: point sources must be re- 
moved, the background moment must be much l ess th an the source moment, and the source 



must fit entirely on the detector. iJeltema et al.l (120051 ) also notes that low S/N images with 



highly symmetric/circular sources will have systematically higher power ratios due to noise. 
A Monte Carlo approach can be employed to estimate the uncertainty associated with noise 
(as described in §3.2). 

In the case of the CLA, the primary limitation is associated with the normalization of 
the two images to have maximum equal intensity. If the images have very disparate S/N, 
the normalization procedure will alter the relative quality of the images. Namely, if the 
maximum intensity of the high S/N image is scaled down to match that of the low S/N 
image, weaker features in the high S/N image would be removed artificially and thus ignored 
in the CLA results. Therefore, the method is best applied to images with similar intensity 
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prior to normalization. As with the PRM, a Monte Carlo approach can be used to measure 
the the effects of noise on the resulting cumulative-distribution function (as described in 
§3.3). 

Unlike the other methods, noise is not a significant source of uncertainty in WTA pro- 
vided that certain criteria are met. We have shown in §2.3.2 that we can limit the affects of 
noise by ignoring pixels with global maxima of a single pixel. Generally, we find that with 
this procedure, the averaged WTA plot is roughly independent of S/N for S/N > 2.0 as long 
as Omax is much greater than 1 pixel (omax > 8 pixels). If emitting regions have S/N and Omax 
below these values, a Monte Carlo approach can be employed to estimate the uncertainty 
introduced by noise. 

For WTA, the principal limitation arises from increased filling factor. In §2.3.2, we 
demonstrated that our ability to identify individual clumps diminishes with greater filling 
factor (see Figure [5]). The reason for this difficulty is that when clumps agglomerate into 
larger structures, it becomes challenging to distinguish the central pixel of that emitting 
region. Thus, isolated clumps in sources with low filling factors are the best candidates for 
WTA to measure indiviudual clump sizes. Large filling factors negatively influence the aver- 
aged WTA results as well: the turnover in the averaged WTA plot becomes less prominent 
with increased emitting surface area in an extended source. Generally, we find that images 
with filling factors above ~30% do not produce a statistically significant amax- 

3. Data Analysis 

In the subsequent sections, we demonstrate the application of the three techniques of 
§2 to a particular SNR, the galactic remnant W49B. 

3.1. Observations and Data Preparation 

We analyzed a 55-ks archival Chandra ACIS observation of W49B (ObsID = 117). The 
remnant was observed with the backside-illuminated S3 chip in the Timed-Exposure Faint 
Mode. Data reduction and analysis was performed using the Chandra Interactive Analysis of 
Observations (ciAo) Version 3.4. We followed the CIAO data preparation thread to reprocess 
the Level 2 X-ray data: we adjusted the charge transfer inefficiency and time-dependent gain, 
and we randomized the PHA distribution and the event positions. ASCA grade 0, 2, 3, 4, 
and 6 events were used in the analysis, and good time intervals were applied. Data below 
0.3 keV and above 8 keV were ignored. Throughout the paper, we assume the distance D 



- 12 - 



to W49B is D = 8 kpc; this value implies an ACTS pixel size of 0.492" ^ 0.019 pc. 

We extracted the global Chandra ACIS X-ray spectrum of W49B (shown in Figure ^ 
using the CIAO command specextract, and ^2.4x10^ total counts were detected from the 
remnant. The X-ray spectrum is quite complex: it has a thermal bremsstrahlung continuum 
and numerous emission features. Prominent blends include both He-like ions and the H- 
like ions of silicon, sulfur, argon, and calcium. Additionally, a strong iron line is detected. 
Modeling of XMM-Newton spectra from W49B suggests at least two plasma components 
are necessary to adequately fit these diverse emis sion features since they are produced at 
different electron temperatures (jMiceli et al. 2006 ). 



Exposure-corrected images of the strong emission features were produced by filtering 
data to the narrow energy ranges of the prominent emission lines (as labeled in Figure [H]): 
Si XIII (1.74-1.93 keV), Si xiv (1.94-2.05 keV), S xv (2.25-2.50 keV), S xvi (2.52-2.71 
keV), Ar xvil (2.98-3.21 keV), Ca xix (3.75-4.0 keV), and Fe xxv (6.4-6.9 keV). We 
multiplied all the exposure-corrected images by the maxima of their exposure maps so that 
every image had units of counts. The spatial distribution of the elements that produce these 
lines is inhomogeneous (Figure (TD] gives examples). The iron is much more prominent in the 
Eastern portion of the remnant, while silicon and sulfur are comparatively more symmetric. 
Additionally, all three ions have discernible X-ray substructure, particularly in the central 
and Western regions where clumping is evident. 

The line emission in these images will have contamination from the underlying thermal 
continuum, and it is necessary to remove this component prior to imaging analyses. We 
estimate the flux of the bremsstrahlung and lines by phenomenologically modeling the X-ray 
spectra from different regions of the remnant. We find that the relative flux of the thermal 
emission in the narrow-band images varies by only ~5% over the entire source. For example, 
bremsstrahlung in the energy range of the Si xill line contributes 35-40% of the flux for all 
regions. Thus, we continuum-subtract the elemental images by reducing the counts in each 
pixel by the average fraction of flux from the thermal component in that spectral band. 



3.2. Ion Symmetry 

We calculate the power ratios P2/P0, P3/P0, and P4/P0 of the continuum-subtracted, 
exposure-corrected images of the five strongest emission features in the W49B spectrum (Si 
XIII, S XV, Ar XVII, Ca xix, and Fe xxv) using the method outlined in §2.1. These five 
images have ~7500-25000 counts each; thus, they satisfy the criteria outlined in §2.4 for 
accurate net moment determinations. The results are listed in Table 1. For each ion image. 
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we determined its centroid and calculated the power ratios using a aperture radius that 
encloses the entire remnant, Rap = 150". To account for the X-ray background, we measured 
the counts per pixel in a circular annulus 155" - 190" from the centroid of each ion image. 
Assuming a constant background, we produced background images for each ion with the 
same dimensions as the source images by setting every pixel to the extracted count-per-pixel 
values. We calculated the moments of the source and of its background separately, and we 
subtracted the background moments from the source moments to obtain the powers using 
these net moments. 

To estimate the uncertainty in th e power ratios , we follow the Monte Ca rlo approach 



employed in the X-ray cluster studies (IBuote fc Tsailll996l : IJeltema et al.ll2005l ). Exposure- 



corrected images (n ormalized to have units of counts) are adaptively-binned using the pro- 



gram AdaptiveBin ( ISanders fc FabianI l200ll ) such that all zero pixels are removed (corre- 
sponding to a minimum of ~2 counts bin~^) to smooth out noise. Then, noise is added 
back in by taking each pixel intensity as the mean of a Poisson distribution and selecting 
randomly a new intensity from that distribution. This process was repeated 100 times for 
each of the five elements, creating 100 mock images per ion. The 90% confidence limits listed 
in Table 1 were obtained using the fifth highest and fifth lowest power ratios from the 100 
mock images of each element. 

We find that the power ratios of Fe xxv are statistically different than the other ele- 
ments. The ratio P2/P0 (a measure of ellipticity) of iron is roughly half the values obtained 
for Si XIII, S XV, Ar xvil, and Ca xix. Similarly, P^/Pq (a measure of ellipticity on smaller 
scales than P2/P0) of Fe xxv is 2-4 times smaller than the lighter ions. However, P3/P0 
(a measure of deviations from mirror symmetry) of Fe XXV is a factor of 4-5 higher than 
the values of the other ions. The surface brightness contours of these ions demonstrate the 
reason for these results (Figure [TTi) . The large P2/P0 and P4/P0 of the lighter ions reflect 
their comparatively elongated distribution relative to Fe xxv. The large P3/P0 of Fe xxv 
arises because of its absence from the Western part of the remnant while it has bright sub- 
structure in the central and Eastern regions of W49B. The other ions have much smaller 
P3/P0 values since their distributions are comparatively more balanced in the Eastern and 
Western portions of the remnant. Therefore, the morphology of Fe xxv is quantitatively 
less elongated and more asymmetric than the lighter ions. 



3.3. Ion Distribution 



We perform the correlation-length analysis (CLA) outlined in §2.2 on images of the five 
strongest emission lines in the Chandra ACIS spectrum. All the ion images have roughly 
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the same maximum pixel intensity (~5 counts) prior to normalization, so the limitations 
of CLA discussed in §2.4 do not apply to the W49B data. We analyzed each combination 
of the five images and produced their correlation-length maps, with examples shown of the 
three most prominent lines (Fe XXV, Si xill, and S xv) in Figure [T21 We find that all ions 
are well-mixed throughout the remnant, except for Fe XX v which is largely absent from the 
plumes of the remnant relative to the other ions. 

Additionally, we produced the correlation-length cumulative-distribution functions shown 
in Figure [121 1-cr confidence ranges of these CDFs were calculated by performing CLA on 
every pair of ten mock images (the same mock images as in §3.2) of each ion. To determine 
whether the ions have similar spatial distributions, we compare the CDFs to those produced 
using auto-correlation (Figure [12]). We find that the silicon-to-sulfur CDF is statistically 
similar to the silicon auto-correlation CDF (with an average difference of O.OOTo") and to 
the sulfur auto-correlation CDF (with an average difference of 0.060cr). This result indicates 
that Si XIII and S xv have the same spatial distribution and are well mixed. Conversely, 
Fe XXV has a very different spatial distribution than Si xill and S XV: the CDFs are very 
disparate from those of the auto-correlations, with average differences of 4.6a from Si xill 
and 3.7a from S xv for correlation lengths under 10". Thus, the iron is more segregated and 
disjoint from the other ions, whereas all except Fe xxv seem to be well-mixed and similarly 
distributed throughout the remnant. 



3.4. X-ray Substructure 

We apply the discrete wavelet-transform analysis outlined in §2.3 for the seven strong 
emission lines in the Chandra ACIS spectrum using Matlab Version 7.4. Figure shows all 
the wavelet-transformed images for five different scales. As per the discussion in §2.2, the 
fields display the X-ray power at the given sizes. At small scales, noise and random fluctua- 
tions dominate, and with increasing Mexican-hat sizes, the distribution and substructure of 
each ion becomes more evident. Consistent with the previous analysis, the Fe xxv emission 
is absent from the Western region of the fields, whereas other ions are more symmetrically 
distributed. 

As in the procedures outlined in §2.4.2, we determine the scale of individual X-ray 
emitting regions as well as the intensity profile in our ACIS ion images. Figure [T^] displays 
examples of identified clumps in the wavelet-transformed images of Fe xxv and their WTA 
plots. The method distinguished successfully the location and scale of the substructure, and 
many different sizes of X-ray emitting regions are evident. Figure [16] shows the averaged 
WTA plot for the five strongest ions in the Chandra ACIS spectrum, Fe xxv. Si xill, S xv. 
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Ar XVII, and Ca xix. We have ignored the global maxima at the smallest scales since these 
are due only to noise and have no physical interpretation. All the ions except Fe xxv have 
the same global maximum, Omax = 17.5" = 0.67 pc. The Fe xxv curve peaks at about 25% 
larger value, with Omax = 22.5" = 0.86 pc and has 40-50% more power than the other ions 
at scales of ~30-60". Thus, the characteristic sizes of X-ray emitting regions of the \ow-Z 
ions are equal, whereas the scale of Fe xxv clumps is greater. Fe xxv contributes much less 
power at small scales relative to these ions as well: below scales of 15", Fe xxv emits 30-60% 
less power than the other ions. These results indicate that the X-ray emitting plasma has a 
characteristic size, and this scale is similar for all ions except for Fe xxv in W49B. We note 
that the ctmax of all ions is sufficiently high to avoid being limited by noise. 

We measure the counts from individual clumps using the CIAO command dmextract 
for a region of area equal to the clump size. We account for the background counts using 
a circular region with diameter 40" located 200" southwest of the remnant's center. We 
calculate the X-ray flux from individual clumps by multiplying the extracted counts by the 
centroid energy of its emission line and dividing by the duration of the observation. In 
Figure [TTl we plot the resulting clump sizes as a function of X-ray flux for Fe xxv. Si xill, 
and S XV. The points at the largest scale reflect the size and flux of these emission lines 
over the entire remnant. On the small scale, the X-ray fluxes go roughly as the cube of 
clump size (i.e., volume). However, when extrapolated to the full size of the remnant, such a 
relationship would severely overpredict the flux of the source. The line fluxes over the entire 
remnant are ^65-100 times less than if they were proportional to volume. The inverse of 
these values gives the fllling factor of the emission in the remnant: 1-1.5%. This small value 
underscores the importance of considering substructure when using X-ray flux to estimate 
remnant properties (such as abundance ratios or elemental masses). 



4. Discerning Heating and Explosion Properties 

From the above analyses, it is obvious that Fe xxv has a distinct X-ray morphology. 
The distribution of Fe xxv is more asymmetric than and disjoint from the other ions, and 
X-ray emitting regions of Fe xxv are larger than those of the lighter ions. Two possible 
scenarios can account for these properties: the iron is insufficiently heated over large areas 
of the remnant, or it was ejected anisotropically during the supernova explosion. From the 
wavelet-transform analysis, we have identified quantitatively the regions with and without 
strong Fe XXV emission, and we can model spectra in these locations to probe plasma 
properties and ultimately, to discern heating versus explosion effects. 
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4.1. Spectral Analysis & Plasma Properties 

We extracted XMM-Newton spectra at seventeen locations identified to have Fe xxv 
clumps (white circles A-Q in Figure [T8l) . W49B was observed three times by XMM: 2004- 
Apr-03 for ~20 ks (ObslD 0084100401); 2004-Apr-05 for ~20 ks (ObslD 0084100501); and 
2004-Apr-13 for 2.5 ks (ObslD 0084100601). We downloaded reprocessed versions of the first 
two observations from the XMM-Newton Science Archiv^. Both observations were taken 
with the PN and both MOS in the Full Window Imaging mode and with the Medium filter. 
After filtering for intervals of high background flaring, the first observation was left with 14.5 
ks of exposure for the PN and 18 ks for each MOS, and the second observation was left with 
14.8 ks for the PN and 18 ks for each MOS. 

For each of the locations circled in Figure [181 we extracted spectra from the PN and 
both MOS detectors using the evselect tool. Our extraction regions were 10" in radius. We 
selected events with PATTERN in the 0-4 range for the PN and 0-12 range for the MOS and 
with FLAG = for both instruments. Response files were generated for each spectrum using 
the rmfgen and arfgen tools to take into account the selection criteria, extraction region size, 
and off-axis angle. Our background spectra were made from a source-free 250" diameter 
circle located ~10' northwest of the center of the remnant. 

Spectra were fit using XSPEC Version 12.4.0. Data were grouped such that a minimum 
of five counts were in each energy bin, and the six MOS and PN spectra from each region 
were fit simulateneously to improve statistics. Standard weighting was used in calculating 
chi-squared values. We modeled the spectra with two components, one cool plasma with 
fixed solar abundances and one hotter plasma with varying supersolar abundances (similar 
to Miceli et al. 2006), in coUisional ionization equilibrium (CIE) using the XSPEC model 
MEKAL (Mewe et al. 1985, 1986; Liedahl et al. 1995). Example spectra and fits from one 
region (circle A in Figure [T8|l are plotted in Figure [191 CIE occurs when the electrons and 
ions reach equipartition, and the electrons are heated only by Coulumb collisions with the 
ions. The timescale r to reach CIE is related to the electron density (in cm~^) and the 
electron temperature Ty (in units of 10^ K; Spitzer 1978): 

T = TrOOn^^Ty^Vears. (12) 

For a young supernova remnant like W49B, with an estimated age of ~1000 years, CIE 
is only realistic with sufficiently high Ue- From radio and infrared studies of W49B, Ug is 
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estimated to be quite high (~1000-8000 cm~^; Keohane et al. 2007), making CIE plasmas a 
plausible scenario. In our fitting of the XMM spectra, CIE models gave chi-squared values 
~10-15% lower than non-equilibrium ionization (NEI) models. All seventeen produced good 
statistical descriptions of the data; Table 2 lists the chi-squared values and degrees of freedom 
associated with each model. 

The best-fit temperatures for the regions with iron emission are given in Table 2: the 
range in the cold plasma temperature is kTi = 0.49-1.19 keV (corresponding to T7 ~ 0.57-1.3^ 
and in the hot plasma temperature is kT2 = 1.81-3.68 keV (corresponding to T7 ~ 2.10-3.96). 



Gener ally, the temperatures increase toward the Eastern side of the remnant. iMiceli et al. 



( 120061 ) found that the strong Si xill line arises from the cold plasma, while the other emission 
lines require the hot plasma. To test this result. Figure [201 shows how the silicon abundance 
(denoted by color) changes as a function of the temperatures of the cool and hot components, 
kTi and kT2. Confidence contours of kTi and kT2 are overlaid to highlight their statistically 
viable values. For constant kTi, the silicon abundance varies less than 10%; by contrast, 
for constant fc T2, the silicon abun dance changes by a factor of :^3-4. Thus, we confirm the 



conclusions of iMiceli et al.l (l2006l ) that the Si xill arises mostly from the cool plasma. 



We extracted XMM-Newton spectra at six locations identified by WTA as weak in Fe 
XXV emission and yet enhanced abundances of lighter ions (the red circles R-W in Figure [T8|) . 
We fit these spectra using a similar model as the regions with strong Fe xxv emission: two 
CIE plasmas, one cool component with fixed solar abundances and one hot component with 
varying supersolar abundances. The resulting chi-squared values and degrees of freedom 
associated with each fit are given in Table 2; all six gave statistically good descriptions of 
the data. To verify the lack of Fe xxv emission in these regions, we measured the reduced 
chi-squared of the fits as a function of iron abundance and compared the results to those 
from regions with strong iron emission. Figure [2T] plots the resulting reduced chi-squared 
values for iron abundances relative to calcium for one region with strong iron emission 
(circle A in Figure [T8l) and for one region with weak iron emission (circle V in Figure [T8|) . 
Since Fe xxv and Ca xix form at the same electron temperatures, we plot the Fe/Ca 
ratio to establish whether iron is comparatively depleted. The iron-weak region has roughly 
constant reduced chi-squared values for iron-to-calcium abundances ranging 10~^ to 0.3 and 
rises sharply at larger ratios. By contrast, fits to the iron-strong emitting regions gave a 
local minimum in reduced chi-squared values at unity. Consequently, the results support the 
WTA identification of these regions as depleted in Fe xxv emission. 

To test whether two plasmas (and not one) are necessary in these areas with depleted 
iron emission, we varied the electron temperatures of the cool component [kTi) and of the hot 
component (/CT2) over a grid of values while letting all other parameters float. This process 
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produced the confidence contours plotted in Figure [221 As the 99% confidence region does 
not intersect the axes (i.e., electron temperatures of zero), plasmas with two temperatures 
are statistically necessary throughout the remnant. 

We model the spectra with weak Fe XX v emission with two CIE plasmas to estimate 
the electron temperatures there. The six regions gave best-fit hot plasma temperatures 
kT2 = 1.80-3.41 keV (corresponding to Tj ~ 2.09-3.96), very close to the values obtained 
for regions with strong Fe XXV emission. At these electron temperatures, Fe xxv should 
have a prominent line (see Figure [25]) if iron is abundant there. However, little-to-no iron is 
X-ray illuminated in these regions, indicating that the iron is depleted from these parts of 
the remnant (with an abundance 1.49±0.88 relative to solar). Thus, insufficient heating is 
not a viable explanation for observed asymmetries, and the anomalous Fe xxv distribution 
must arise from an anistropic ejection of iron. 



4.2. Explosion Mechanism 



Since we have determined that the asymmetries in W49B result from an anisotropic 
ejection, we now explore supernova explosions that could lead to the observed Fe xxv mor- 
phology. Bipolar explosions, a special class of core-collapse supernovae, are highly asymmet- 
ric explosions that would distribute heavy metals similarly to those in W49B. Core-collapse 
supernovae occur when the iron cores of massive stars (M ~ 8 — I3OM0) collapse to form 
neutron stars or black holes. These explosions are categorized based on whether the pro- 
genitor lost its hydrogen envelope (Type lb or Type Ic) or not (Type II). A core-collapse 
supernova could produce a bipolar explosion if the progenitor's core is undergoing sufficiently 
large rotatio n. Numerous mechanisms could be responsible for the resulting bipolar nature 



of the event ( Ramirez- Ruiz et al 



2002 



MacFadven fc Wooslevl[l999[: [MacFadven et al. 



2001 



Duncan fc Thompson 



Komissarov fc Barkov 



992I: [Wheeler et al.[[2000[ : [Bucciantini et al.[[2007[ : [Burrows et al.[[2007[ : 



2003) 



Typically, core-collapse SNe have kinetic energies around 10 ergs. However, in the last 
ten years, numerous detections of hypernovae — supernovae with high ejecta velocities and 
sometimes large kinetic energies (e.g., Nomoto et al 2001) — have challenged this paradigm. 
Some of these hypernovae may originate from bipolar explosions (Maeda & Nomoto 2003). 
With the number of bipolar explosion/hypernovae candidates growing, the distinguishing 
properties of these extreme sources are being revealed. In the case of SN 1998bw, the [Fe 
11] lines were broader than the [O l] AA doublet, indicating that th e ^^Ni (which decay s 
into iron) was ejected at larger velocities with respect to the oxygen (iMazzali et al.|[200l[ ). 
A spherically symmetric model of the collapse could not account for this observation, and 
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Maeda et al.l (|2002| ) successfully reproduced the result using an axisjTiimetric model. In this 



scenario, the kinetic energy is assumed to decrease as a function of solid angle from the 
polar axis, causing the nickel to be more efficiently synthesized and ejected at high velocity 
in that direction. Perpendicular to the rotational axis, mechanical heating decreases; thus, 
this phenomenon predicts that lower-Z elements are ejected more isotropically and at smaller 
velocities than the nickel (e.g., Mazzali et al. 2005). 

Generally, these asymmetries alter the chemical yield of bipolar explosions relative to 
typical CC SNe. While the typical core-collapse event ejects ~ 0.07 — O.ISM© of ^^Ni, bipolar 
explosions can have much gr eater nickel yields that als o increase with asphericity, progenitor 
mass and explosion energy (jUmeda fc Nomotd l2008l ) . Examples of nickel masses derived 
for bipolar explosion candidates include ~0. 25-0. 45 Mq (2003dh: Mazzali et al. 2003), 
~0.45-0.65 Mq (20031w: Mazzah et al. 2006), and ~0.2-0.7 Mq (1998bw: Woosley et al. 
1999; Kaneko et al. 2007). 

The abundance ratios obtained in our spectral fits can be u sed to const a in wh ether 
W49B may be the remnant of a core-collapse, bipolar explosion. iMiceli et al.l (120061 ) per- 
formed similar analyses for the bright, central region of W49B (circle I in Fig. [T5]l and 
cor npared its abundance rat ios to aspherical explosion models of extremely metal-poor stars 
by iMaeda fc Nomotd (120031 ). Based on low silicon and sulfur abundances relative to iron, 
these authors exclude massive progenitors (Mzams > 25 Mq) as well as large explosion ener- 
gies {E > 10^^ erg). This conclusion assumes that the abundances in the central region are 
refiective of the entire remnant. However, the iron abundance is greatly enhanced in that 
location: as listed in Table 2, this area (circle I) has the lowest Si/Fe and S/Fe abundance 
ratios of all twenty-three regions in Figure [181 Additionally, we have demonstrated in §3 that 
the iron has statistically different symmetry, distribution, and substructure than the other 
elements. This attribute makes it imperative to average over many regions of the remnant 
to get an accurate representation of the abundance ratios. However, it is not adequate to 
estimate abundances from fits to the total spectrum because the central bar dominates the 
emission, giving uncharacteristicly low values for the lighter elements relative to iron. For 
example, the total spectrum gives a Si/Fe ratio of ~0.18 and a S/Fe ratio is ~0.07, close 
to the circle I values. We also note that the Miceli ratios are derived solely from the hot 
plasma, which affects the accuracy of the silicon abundance determination since it comes 
mostly from the cool plasma (as shown in Figure [20]) . Consequently, the Miceli estimate of 
the silicon abundance is low since it does not include the primary heating source for that 
element. 



To obtain an accurate measurement of the global abundance ratios in W49B, we fit the 
XMM spectra from the twenty-three regions in Figure [181 seventeen with strong Fe xxv 
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emission and six with depleted Fe xxv emission. We plot the mean abundance ratios of the 
hot plasma (since the cool plasma is fixed to solar values) relative to iron for silicon, sulfur, 
argon, and calcium in Figure [21] as well as list their values and 1-a errors in Table 2. The 
ratios have fairly large dispersions, with all spanning factors of 5-8. This result demonstrates 
that relative abundances of the elements change throughout the remnant, underscoring the 
importance of averaging over many regions in this calculation. We increase the silicon abun- 
dances plotted in Figure [2^ by 30% from their Table 2 values to account for the contribution 
of the cool component. We estimated this factor by re- fitting the regions' spectra and fix- 
ing all elements of the cool plasma to solar metallicity except the silicon abundance, which 
was allowed to vary. Following this procedure, the combined silicon abundance of the two 
plasmas was 30% greater than the hot component alone. 

In Figure we compare our abundance ratio s to those predicted by six core- collapse. 



Type Ic models: four energetic, spherical explosions (INomoto et al.ll2006l ) and two aspherical, 
bipolar explosions (25A, 25B: Maeda & Nomoto 2003). The aspherical cases have the best 
agreement with our measured abundance ratios; both models' values are within the 1-a 
range for all elements in the W49B spectrum. The primary difference between the two 
aspherical models is their jet opening angles 6: 25A has 6 = 15° and 25B has 6 = 45°. For 
the fixed kinetic energy, a larger 6 means less mechanical energy is released per unit solid 
angle in 25B, causing the explosion to appear less luminous and altering the nucleosynthetic 
products. Despite these distinctions, the dispersion in measured abundance ratios for W49B 
prevents us from discriminating between these two alternatives. Nonetheless, we can exclude 
the spherical explosion cases based on the silicon and sulfur ratios. Therefore, we find that 
a bipolar origin with typical or slightly larger kinetic energies of W49B is consistent with its 
global abundance ratios. 



4.3. Mass Estimates 

The explosion mechanism can also be constrained by estimating the mass of metals 
in W49B. We perform this calculation using the elements' filling fa ctor (as determined by 



the WTA) and the strength of the emission lines. Appendix B of iLazendic et al.l (120061 ) 
summarizes a similar method to find masses of individual regions of a source; we adapt this 
technique to the full W49B remnant and account for substructure using WTA. 

Masses are determined based on parameters given in the X-ray spectral fits. For a given 
electron temperature, the fiux of an emission line is proportional to its emission measure, 
£^Miine = / ngriidV, where rij is the ion number density and dV is the emitting volume. If ( 
is defined as the ratio of electron and ion number densities, Ue/ni, then rii is given by 
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^.= (^j , (13) 

where we have set dV = V, the volume of a chosen emitting region. From there, it is simple 
to derive the ion mass M, over the remnant: 

M, = fn^Vm, = fm, (^^^^ (14) 

where vni is the mass per ion in grams and / is the filling factor of the ion's emission over 
the remnant. 

The primary uncertainty in the above mass calculation is the electron number density 
Tif.. However, since the plasma is in a collisional equilibrium state, we can set a physi- 
cally motivated lower limit on ng. By setting r to the age of the remnant in Equation Wl\ 
(~ 1000 years) and setting the electron temperature to the greatest value obtained by our 
spectral fits, T7 ~ 3.96, we find ne > 60 particles cm~^. Additionally, two reasonable sce- 
narios define a strict range for C, = rie/rii (see §3 of Hughes et al. 2003 for a discussion). 
In the extreme case of a pure metal plasma, there are 50 electrons per Fe ion {( = 50). 
Conversely, if the plasma has solar abundance (a large hydrogen content), there are 3x10^ 
electrons per Fe ion {( = 3x10^). An intermediate case, one where hydrogen mass is equal 
to that of metals (where ( = 600), is also plausible. Our uncertain ties are defined by th e 



range of possible Ue, and our values are consistent with the results of iKeohane et al.l (120071 ). 



We can determine the filling factor / using our WTA results and direct observables. 
The filling factor is defined as the ratio of the summed volume of individual clumps 
to the total volume of the entire source, Vt- Assuming the distance D = 8 kpc, we esti- 
mate the radius of X-ray emission in W49B is i? ~ 138" = 5.3 pc = 1.65 xlO^^ cm, giving 
Vt = |7ri?3 = 1.9 X lO^s cm . We find the summed volume of individual clumps by multiply- 
ing the mean volume of clumps from our WTA, Knean = f^ra^g^, by the number of clumps, 
Nc- The number of clumps is determined simply by computing how many clumps with Vmean 
and a fiux Fc are necessary to account for the total fiux Ft of the source (as calculated using 
the procedure in §3.3: 

Nc=^. (15) 



We opt to measure Nc this way instead of counting all the clumps in the WTA images 
because the latter would require us to define aw/a threshold for an individual clump. 
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Putting all of these relations together, the filling factor / is 



Thus, 



Vc A^cKiean -^T ^max ( a\ 



Ft 3 /47rEMune\i/2 



M.^^aLx-.(^^) . (17) 



We use Equation [17] to calculate the total mass of iron in W49B. From Figure 
Omax ~ 22.5" = 0.86 pc for iron, and from Figure fTTl = Ft/Fc = 3.75 for iron. These 
values give a filling factor / = 0.016, an order of magnitude smaller than previo us estimates 



that were based on assumptions about the extension of the emitting plasma (IMiceli et al. 



20061 ) . Thus, a rigorous measurement of / demonstrates that the value is actually much 



smaller, and iron especially must emit its power in concentrated locations of the remnant. 

We find the emission measure of iron EMij-o^ by the following procedure. Using the 
model from §4.1, we fit the global Chandra ACIS X-raj spectrum of W49B (Figure [9]) to 
obtain the emission measure of the hot-plasma component: EM^ot = 1-805 x 10^^ cm~^. 
To compute the emission measure of an individual line EMu^e, we reduce EM^^ot by the 
ratio of fiux in the line to that of the entire spectrum. X-ray fiuxes are computed using the 
XSPEC command flux over the energy ranges of interest for a phenomenological model of 
the lines and continuum. Based on this procedure, iron contributes 23% of the fiux to the hot 
component, so EMi^on = 4.29 x 10^^ cm~^. 

Using the equations above, we find an iron mass of Afiron ~ 0.12 Mq assuming a solar 
abundance plasma. For the case of hydrogen mass equal to that of metals, we obtain an iron 
mass of Miron ~ 0.80 Mq. The condition of a plasma in coUisional equilibrium defines the 
upper-limit iron mass in W49B (since it gives the minimum Ue necessary to reach CIE by the 
age of the remnant): Miron ~ 1-29 Mq. While this range of iron masses is large (refiecting 
the uncertainty in electron density), masses above 1.29 Mq and below 0.12 Mq can be 
excluded because it would require physically implausible metal abundances. Furthermore, 
the large iron mass cannot be attributed to a thermonuclear explosion since the abundance 
ratios relative to iron in Figure [2l] are much greater than those expected from Type la 
events. Although the derived iron masses are larger than those usually predicted for core- 
collapse supernovae, the lower limit is close to those values for typical bipolar explosions 
(^ O.I5M0). Additionally, as noted in §4.2, recent models of bipolar explosions demonstrate 
that nickel yields (which decays into iron) can actually be much greater than the typical core- 
collapse scenario (up to 4 Mq) with increased asphericity, explosion energy, or progenitor 
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mass (Umeda & Nomoto 2008). The anisotropic ejection of iron (as shown in §4.1) and 
abundance ratios (given in §4.2) in W49B are consistent with such an asymmetric explosion 
that would overproduce nickel. Generally, we emphasize that this calculation is an order-of- 
magnitude estimate, and precise measurement of the filling factor is an important step in 
obtaining a plausible range of masses. 

A few factors contribute to the uncertainty in our analyses. As noted previously, the 
largest source of error is the electron number density since we do not know precisely the 
plasma's composition. Wc attempt to limit this uncertainty by defining a plausible range in 
its value. Another source of error is the physical modeling of the X-ray spectrum. While the 
CIE model produces a chi-squared value lower than an NEI plasma, it is barely a statistically 
significant difference (~10%). However, NEI models give abundance ratios within the errors 
of our CIE results, so we expect this uncertainty is minimal. Other sources of error are 
instrumental in nature: the off-axis point-spread function (PSF) effects with Chandra and 
the possible clump contamination with XMM-Newton from its limited spatial resolution. 
The PSF of Chandra ACIS increases as a function of energy and of off-axis angle. At its 
aimpoint, the PSF is ~1" at 1.49 keV and ~2" at 6.4 keV; at a distance 6' off-axis, the PSF 
increases to ~3" at 1.49 keV and ~5" at 6.4 keV (where we have characterized the PSF as 
the scale where the encircled energy fraction is 90%). We anticipate that these effects do 
not alter our results much since the overall extent of W49B is small (~4' in diameter) and 
the characteristic clump sizes of all ions are many factors greater than these values. We 
minimize the uncertainty associated with XMM-Newton clump contamination by selecting 
emitting regions that were sufficiently isolated from other substructure. 

The SNR mass estimate in this paper is the first with the contribution of substructure 
rigorously defined. Measurement of ejecta masses requires detailed knowledge about the scale 
of the emitting regions since emission measure increases linearly with volume. Previous SNR 
studies assume generally that a source's emission is homogeneous (e.g., Hughes & Singh 1994) 
or estimate the filling factor from basic arguments about the plasma's extension (e.g., Miceli 
et al. 2006) when calculating physical parameters. However, the complex ion distribution 
and low filling factor in W49B demonstrates the necessity for more rigorous approaches. 
Consequently, it is vital to probe accurately the role of substructure when determining 
global properties of a supernova remnant. 

5. Summciry 

In this paper, we have explored three methods that quantify the morphologies of ex- 
tended sources. We have demonstrated the accuracy of these techniques on relevant synthetic 
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data, and we have exemplified the utihty of these tools by applying them to the Chandra 
ACIS observation of the SNR W49B. 

Using a multipole power-ratio technique, we measured the asymmetry of the ion distri- 
butions in W49B. We discovered that the morphology of Fe xxv is less elongated and more 
asymmetric than the lighter ions. We applied an adapted two-point correlation method 
called correlation-length analysis to compare the ions' spatial distributions. We determined 

that the Fe XXV emission is more segregated and disjoint from that of the lower-Z ions, 
which are well-mixed and similarly distributed throughout the remnant. Additionally, we 
measured properties of the X-ray ions' substructure in W49B using wavelet-transform anal- 
ysis. We found that the mean clump scale of Fe xxv is roughly 25% larger than those of the 
hghter ions. Together, these results demonstrate quantitatively that Fe xxv has a distinct 
X-ray morphology. 

Two potential scenarios have been proposed to explain the anomalous properties of 
iron in W49B (Miceh et al. 2006; 2008; Keohane et al. 2007): insufficient heating of 
iron over large areas of the remnant or anisotropic ejection of iron during the supernova 
explosion. By fitting XMM-Newton spectra in twenty-three locations with strong and weak 
Fe xxv emission (as identified by our wavelet-transform analysis), we determined that the 
electron temperatures throughout the remnant should heat any iron there to radiate at X-ray 
temperatures. Thus, iron is strongly depleted from the identified regions, and the distinct 
X-ray morphology must arise from an anisotropic ejection. Therefore, any models Tised to 
describe the explosion mechanism of W49B must include a means for the iron to be ejected 
an asymmetric way. 

We further constrained the nature of the explosion by measuring the mean abundance 
ratios of elements in our twenty-three X-ray spectral fits. The ratios varied tremendously in 
different locations, demonstrating the importance of averaging the ratios over many regions 
to obtain the global abundances in W49B. We found that the mean ratios are consistent 
with bipolar explosion models of massive stars. This type of explosion could also account 
for the anistropic ejection of nickel since heavy elements are preferentially ejected along the 
polar axis of the progenitor in this scenario. We calculated the total iron mass in W49B 
by determining the filling factor (using our wavelet-transform analysis) and the emission 
measure. We found that the iron mass and abundance ratios are broadly consistent with a 
core-collapse supernova. 

This work is a first attempt to rigorously describe the physical properties of supernova 
remnants. SNRs play a vital role in our galaxy, and the wealth of data available on these 
sources is tremendous. Robust methods for comparison within and between these complex 
targets can advance SNR science as well as probe many poorly understood phenomena. We 
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plan to apply the techniques presented here to other sources to constrain further processes 
as particle acceleration, interactions with environments, and explosion mechanisms. 

We thank K. C. Schlaufman and S. E. Woosley for helpful discussions and contributions. 
This work is supported by DOE SciDAC DE-FC02-01ER41176 (LAL and ER-R) and a 
National Science Foundation Graduate Research Fellowship (LAL). 
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Table 1. Power Ratios and Associated 90% Confidence Limits 





P2/P0 




Ion 


(x 10-^) 


(x 10-7) (x 10-7) 



Fexxv 72.68116-95 O.SSto^s 

Si XIII 137.54tJ7.47 2.87tJ-^« 0.92t°-^8 

Sxv 141.64+1^:94 4_i5+i:o9 145+0:55 

Ar XVII 141. 29t?^-^^ 3.45^}^^ l-Slto rg 

Caxix 136.19t?7:o7 3.42t}J2 i.Tet^J^ 



Table 2. Best-Fit Temperatures and Hot-Plasma Abundance Ratios (Relative to Iron by 

Mass) 



Region fcTj kT-z Si/Fe S/Fe Ar/Fe Ca/Fe x^/d.o.f. 

(keV) (keV) 



Regions with Strong Iron Emission 



A 


0.91±0.10 


2.92±0.31 


0.46±0.14 


0.26±0.08 


0.08±0.03 


0.08±0.03 


939.7 / 918 


B 


0.59±0.11 


1.99±0.12 


0.31±0.06 


0.19±0.03 


0.05±0.01 


0.05±0.01 


482.5 / 514 


C 


1.05±0.11 


3.33±0.33 


0.57±0.22 


0.56±0.22 


0.12±0.05 


0.14±0.06 


694.3 / 701 


D 


0.70±0.10 


2.51it0.17 


0.18±0.06 


0.12±0.03 


0.04±0.01 


0.03±0.01 


877.9 / 934 


E 


0.67±0.17 


2.43±0.17 


0.40±0.10 


0.32±0.07 


0.08±0.02 


0.07±0.02 


832.7 / 855 


F 


1.19±0.15 


3.68±0.38 


0.33±0.21 


0.17±0.08 


0.09±0.04 


0.08±0.04 


461.5 / 480 


G 


0.49±0.14 


2.46±0.23 


0.27±0.09 


0.14±0.04 


0.05±0.02 


0.04±0.01 


387.9 / 416 


H 


0.79±0.13 


2.55±0.13 


0.30±0.06 


0.23±0.04 


0.08±0.02 


0.06±0.01 


1358 / 1405 


I 


0.80±0.13 


2.54±0.13 


0.13±0.04 


0.14±0.02 


0.04±0.01 


0.04±0.01 


1366 / 1410 


J 


0.66±0.29 


1.81±0.32 


0.28±0.17 


0.15±0.07 


0.08±0.04 


0.06±0.03 


351.3 / 388 


K 


0.60±0.34 


1.86±0.12 


0.39±0.08 


0.30±0.05 


0.07±0.02 


0.07±0.02 


663.9 / 676 


L 


0.95±0.11 


3.02±0.34 


0.33±0.14 


0.50±0.21 


0.25±0.11 


0.21±0.09 


617.6 / 657 


M 


0.97±0.11 


2.18±0.29 


0.77±0.46 


0.51±0.31 


0.17±0.10 


0.14±0.08 


644.7 / 738 


N 


0.91±0.21 


2.30±0.22 


0.39±0.10 


0.29±0.09 


0.07±0.03 


0.06±0.02 


819.3 / 868 





0.92±0.28 


1.93±0.18 


0.44±0.12 


0.23±0.07 


0.07±0.02 


0.06±0.02 


863.9 / 941 


P 


l.OOitO.13 


2.42±0.32 


0.70±0.33 


0.50±0.26 


O.llitO.06 


0.12±0.06 


820.0 / 818 


Q 


0.69±0.11 


2.28±0.53 


0.31±0.13 


0.23±0.10 


0.03±0.01 


0.08±0.03 


207.9 / 245 


Regions with WeaJf Iron Emission 


R 


0.79±0.11 


3.41±1.03 


0.56±0.23 


0.47±0.19 


0.38±0.15 


0.22±0.09 


217.0 / 271 


S 


0.89±0.25 


2.17±0.45 


0.53±0.37 


0.37±0.25 


0.07±0.05 


0.08±0.06 


302.0 / 339 


T 


0.50±0.16 


1.80±0.19 


0.55±0.24 


0.53±0.17 


0.12±0.05 


0.13±0.04 


327.7 / 376 


U 


1.02±0.07 


2.40±0.37 


0.26±0.23 


0.49±0.44 


0.23±0.20 


0.26±0.23 


558.1 / 627 


V 


0.55±0.14 


1.85it0.23 


0.65±0.28 


0.40±0.15 


0.11±0.05 


0.08±0.04 


276.9 / 319 


W 


0.83±0.15 


2.16±1.11 


0.70±0.52 


0.42±0.31 


0.09±0.07 


0.06±0.04 


242.5 / 283 
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Fig. 1. — A composite image of W49B as observed with Chandra ACIS in the soft band 
(0.3-1.50 keV; blue), medium band (1.5-2.5 keV; green), and hard band (2.5-8.0 keV; red). 
The image has been smoothed using a Gaussian kernel of width cr = 3 pixels along each 
axis; the scale bar is 100 pixels long. The remnant has a bright central bar, with prominent 
plumes on both sides. The hard X-ray emission generally appears enclosed by the soft- 
and medium-band emission, except for the prominent hard X-ray emission in the Northwest 
region of the remnant. Soft-band emission is weak compared to the other bands because it 
is prefentially absorbed through the Galactic plane. 
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Group 1: 

50 clumps with 

Fixed Position & 

Sizes Selected from 
a Gaussian Dist. 



Group 2: 

50 clumps with 

Fixed Position & 

Larger Variation 
in Size 



Group 3: 

50 clumps with 
Random Position & 
Smaller Variation 
in Size 



•I* i I •• 



Fig. 2. — Examples of synthetic data used in CLA. The first set of ten synthetic images 
(Group 1; examples in top row) contains 50 clumps of fixed position with sizes chosen from a 
Gaussian distribution with mean /igizc = 10 pixels and a standard deviation o"sizo = 3 pixels. 
The second set of ten images (Group 2; examples in middle row) is identical to the first group 
except the standard deviation in size was increased (to agizc = 4 pixels). The third set of ten 
images (Group 3; examples shown in bottom row) has 50 clumps with randomly-assigned 
positions and sizes chosen from a Gaussian distribution with smaller standard deviation, 
(Tsize = 1 pixel. The scale bar is fifty pixels in size. We performed CLA on every combination 
(ninety different pairings) of the ten images in each group and produced the three mean 
correlation-length CDFs and their 1-a ranges shown in Figure [31 
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Fig. 3. — Correlation-length cumulative-distribution functions (CDFs) and their 1-cr confi- 
dence ranges for the three groups of synthetic data (Group 1: blue; Group 2: red; Group 3: 
green). Since the CDFs of Groups 1 and 2 overlap, the two have statistically similar spatial 
distributions, indicating that CLA is not sensitive to moderate changes in clump size. By 
constrast. Group 3 has a very distinct CDF curve since its clumps have randomly-selected 
locations (unhke Groups 1 and 2). These results suggest that CLA is a good means to de- 
termine if images have similar distributions, regardless of whether the clumps are the same 
size. 
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Fig. 4. — An example of wavelet-transform analysis, a: Synthetic data of a single two- 
dimensional Gaussian function with a width a = 10 pixels, b: Plot of the relative power 
(w/a) as a function of Mexican-hat radius a for the central pixel of the clump. Increased 
values of w/a indicating that more power is emitted at a given scales. Thus, local and 
absolute maxima denote the size of an emitting region. Indeed, the peak of the plot occurs 

at flmax ~ ^- 
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Fig. 5. — Synthetic datasets (left) with different numbers of randomly distributed Gaus- 
sians of equal width a = 10 pixels: = 30 clumps, 100 clumps, and 300 clumps. All three 
arrays are the same size, 300 x 300 pixels, and the scale bar is fifty pixels in size. An increase 
in corresponds to emission over more surface area (increased filling factor), so the power 
at larger scales will rise as a function of A^. The averaged WTA plot ( {w)/a versus a; right, 
solid lines) reflects this relationship: the A^ = 30, 100, and 300 images had Omax of 16 pixels, 
22 pixels, and 26 pixels, respectively. Additionally, the maxima become less prominent with 
larger A^ because of this increased filling factor. Thus, the peak of the aveaged WTA plot 
is sensitive to the filling factor of an emitting region. The averaged WTA plot for isolated 
clumps (right, dashed lines) accurately reflects the scales of those structures. The fraction 
of identified isolated clumps decreases with larger filling factors: 24 isolated clumps were 
identified in the A^ = 30 clumps case, whereas 43 isolated clumps were identified in the A^ 
= 300 case. 



-34- 





1 2 4 6 8 10 15 

a (pixels) 



Fig. 6. — Synthetic data of pure noise (left) and the resulting averaged WTA plot (right). 
The scale bar on the synthetic data is fifty pixels in size. The intensity of each pixel with 
noise was chosen from a Poisson distribution with standard deviation o"noise = 3 counts. The 
power profile peaks at Omax = 1 pixel, indicating that noise contributes the most power at 
the single pixel scale. 
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Fig. 7. — Examples of synthetic data of pure noise and 50 clumps combined (a and b) and 
their WTA results (c). The clump sizes were chosen from Gaussian distributions with a mean 
/^sizes = 30 pixels and standard deviation cxsizes = 6 pixels (a) or o"sizes = 10 pixels (b). The 
averaged WTA plot (c) for the datasets overall (shown with solid lines) have global maxima 
at ttmax = 1 pixel because of the noise. However, if we ignore pixels with global maxima 
at cimax = 1 or 2 pixels, the contribution from noise is drastically reduced and the global 
maxima become cimax/Usizcs = 30 pixels. The method extracted successfully the characteristic 
size of the clumps, regardless of their variance in scale. 
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Fig. 8. — Synthetic data of noise and clumps combined (Figure 9a), convolved with Gaussian 
functions of widths a = 2, 4, 6, and 20 pixels to simulate PSF effects on the clumps in an 
astronomical image. Generally, the scale of the convolving Gaussians does not shift the WTA 
peak, and larger Gaussians decrease the power contribution of noise at the smallest scales. 
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Energy (keV) 

Fig. 9. — Background-subtracted global X-ray spectrum of W49B obtained with Chandra 
ACTS. The spectrum is composed of a thermal bremsstrahlung continuum and many promi- 
nent emission lines, including blends of Si, S, Ar, Ca, and Fe. No emission features are 
detected below ^1.7 keV because the source is heavily absorbed since it is viewed through 
the Galactic plane. 



Fe XXV 



Si XIII 




Fig. 10.— ACIS X-ray images of Fe XXV (left). Si XIII (middle), and S XV (right). Images 
are ^250" across and were smoothed using a Gaussian kernel of width a = 3 pixels along 
each axis. Iron has a very distinct spatial distribution compared to the lighter ions, and all 
have many X-ray substructures. 
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Fig. 11. — Surface brightness contours overlayed on smoothed images of Fexxv (top), Si 
XIII (middle), and S XV (bottom). Contours mark 25% (yellow), 50% (orange), and 75% 
(red) of the maximum surface brightness in each image. The green circles mark the centroids 
used in the power-ratio analysis; the centroid of iron is offset from the other elements by 
~25". The surface brightness contours reflect the lack of iron in the Western portion of 
the remnant, while the other elements are comparatively more symmetrically and diffusely 
distributed. The power ratio method quantifies these morphologies and statistically confirms 
the asymmetry of iron relative to silicon, sulfur, argon, and calcium. 
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Fig. 12. — Example maps of CLA results: silicon versus iron (a); sulfur versus iron (b); 
silicon versus sulfur (c). Maps are 300" across. Yellow and red signify large correlation 
lengths, i.e. regions in which the two given ions do not have statistically similar spatial 
distributions. Iron appears to be very segregated from silicon and sulfur, while the latter 
seem to be well-mixed, a result that is confirmed quantitatively in Figure [ 
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Fig. 13. — CLA results for iron, silicon, and sulfur. Solid lines and their respective shaded 
regions denote the cumulative distribution functions and 1-a confidence ranges for iron auto- 
correlation (blue), silicon auto-correlation (red), and sulfur auto-correlation (green). Dashed 
lines and their respective shaded regions indicate the CDFs and 1-cr ranges for silicon vs. 
iron (pink), sulfur vs. iron (light blue), and silicon vs. sulfur (yellow). The silicon versus 
sulfur case is similar to those of the ions against themselves, indicating the two elements have 
statistically similar spatial distributions. By contrast, the CDFs of those elements against 
iron are very disparate from the other curves, demonstrating quantitatively the segregation 
of iron relative to the other ions. 
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Fig. 14. — Wavelet-transformed images of the seven ions detected in the Chandra ACIS 
observation of W49B for four different sized Mexican hats (a = 1", 2", 8", 15", and 30"; 
recall 1 pixel = 0.492"). Images are normalized within each column, and the scale bar is 
50" ~ 1.93 pc in size, for a distance D = 8 kpc. Red and yellow signifies a lot of emission at 
those locations and scale, while blue denotes little-to-no power at those locations and scale. 
Raw images are continuum subtracted prior to the WTA (see text for details). 
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Fig. 15. — Above: Two wavelet-transformed images (w/a) of Fe XXV for Mexican-hat sizes 
a = 5"and 13". The scale bar is 50" ~ 1.93 pc in size. White boxes identify the clumps with 
WTA plots underneath. Below: Example WTA plots {w/a versus a) for individual clumps 
identified at each scale. The x-axis is given in units of arcseconds and parsecs. Errors are 
calculated via WTA assuming Poisson statistics. As in the synthetic data analysis, absolute 
maxima occur at the scale of most power, i.e., the size of an emitting region. Many Fe xxv 
clump sizes are evident. 
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Fig. 16. — The averaged WTA plot for the five strongest ions in the Chandra ACIS spectrum 
of W49B. Each point represents the power measured at the given scale for Fe xxv, Si xill, 
S XV, Ar XVII, and Ca xix. Curves are normalized to unity for easy visual comparison, and 
we ignored pixels with absolute maxima at the smallest scales to avoid contribution from 
noise. We estimate errors to be negligible (~ 10~^) since each point represents the average 
value of ~150,000 pixels. The absolute maxima are at amax = 22.5" ~ 0.86 pc for iron and 
c^max = 17.5" ~ 0.67 pc for all other elements (where 0.492" = 1 pixel). Iron also contributes 
40-60% less power at scales <15" compared to the other ions. The turnover at very high a 
values occurs as the Mexican-hat scale approaches the overall size of the remnant. We note 
that the Omax for all elements is much greater than the on- and off-axis point-spread function 
of ACIS, so we are not limited by instrumental resolution. 
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Fig. 17. — X-ray flux versus clump size for Fe xxv (red circles), Si xill (blue triangles), and 
S XV (green stars). The points at the largest scale represent the flux and size of the entire 
remnant. The orange line represents the case if flux went as the cube of clump size (i.e., 
volume) . 
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Fig. 18. — Smoothed image of the iron emission in W49B. White circles identify regions with 
iron clumps (as identified by the wavelet-transform analysis) where we extracted XMM- 
Newton spectra. Red circles denote regions with weak iron emission where we extracted 
XMM-Newton spectra. Table 2 lists the best-fit electron temperatures and abundance ratios 
by mass relative to iron for each region. 
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Fig. 19. — Example XMM-Newton spectra, fits, and residuals for circle A in Figure 18. 
Best-fit electron temperatures and abundance ratios by mass relative to iron are given in 
Table 2. By fitting the six XMM spectra from circle A simultaneously (each plotted in 
a different color), we improve our count statistics dramatically, and all spectra fits have 
reduced chi-squared values within 10% of unity. For the model in this figure, x^/d-o.f. = 
1366/1410. 
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Fig. 20. — Dependence of the silicon abundance on the temperature of the cool plasma {kTi) 
and hot-plasma {kT2) relative to solar. For constant kTi, the abundance only changes by 
~10%, while for constant kT2, the abundance changes by a factor of 3. Based on these 
resulst, we conclude that the silicon emission arises mostly from the cool plasma. 
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Fig. 21. — Plot of reduced chi-squared as a function of iron abundance relative to calcium 
(with respect to solar by number). While the region with strong iron emission (circle A) has 
an absolute maximum at unity, the region with weak iron emission (circle V) has roughly 
constant chi-squared values for five orders then a sharp incline. This result statistically 
confirms the iron-depleted locations identified by WTA. 
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Fig. 22. — 68%, 95%, and 99% confidence contours for kTi and kT2 in a region with weak 
iron emission (as identified by our WTA analysis; circle V in Figure [T8l) . As the 99% 
confidence region does not intersect with either axis, two plasmas with different temperatures 
are required throughout the remnant. 
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Fig. 23. — Plot of relative strength of Fe ions as a function of electron temperature Te in 
a CIE plasma. The best-fit models of the regions with weak iron emission give electron 
temperatures corresponding to those where the Fe xxv should have a prominent emission 
line. However, there is little-to-no Fe xxv emission in these regions. As the plasma is 
sufficiently heated to irradiate iron, this result indicates that iron is greatly depleted in these 
regions. Thus, the anomalous distribution of iron must arise from an anisotropic ejection of 
nickel during the supernova explosion. 
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Fig. 24. — Plot of mean abundances (by mass) and their range of values relative to iron 
of silicon, sulfur, argon and calcium. See Table 2 for the abundance ratios relative to iron 
for all twenty-three regions. For comparison, we show the abundance ratios predicted by 
models of four spherical explosions (from Nomoto et al. 2006) and two aspherical explosions 
(from Maeda & Nomoto 2003) as well. The dispersion in abundance ratios, particularly 
silicon and sulfur, demonstrate the importance of averaging abundance ratios over the entire 
remnant. Generally, the spherical models are not consistent with the measured ratios, while 
the aspherical models succeed at fitting our obtained values, evidence that W49B is probably 
from a bipolar explosion. 



